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We have investigated the properties of the resistive state of the narrow superconducting channel 
of the length L/<f = 10.88 on the basis of the time-dependent Ginzburg-Landau model. We have 
demonstrated that the bifurcation points of the time-dependent Ginzburg-Landau equations cause 
a number of singularities of the current- volt age characteristic of the channel. We have analytically 
estimated the averaged voltage and the period of the oscillating solution for the relatively small 
currents. We have also found the range of currents where the system possesses the chaotic behavior. 
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It is well-known that the resistive state appears in a su- 
perconductor in the range of currents jc < j < J2 • When 
starting from the superconducting state and increasing 
the current, it is obvious that the superconducting state 
switches to the resistive one at a certain critical current. 
Further increase of the current leads to the switching 
from the resistive to the normal state at the upper criti- 
cal current J2. In the resistive state the superconductivity 
and a static electric field coexist in a system. It causes 
the appearance of phase slip centers (PSCs). The phase 
of the order parameter (OP) in the resistive state period- 
ically drops by 27r in a set of sample points. This process 
leads to the oscillation of the amplitude of the OP. It has 
been shown that PSCs generate certain singularities on 
the current- voltage characteristic (CVC). It possesses a 
stair like structure [l|, |2| . 

The properties of the resistive state of different types 
of the superconductors have been investigated recently. 
The influence of the ratio between the relaxation times of 
the time-dependent Ginzburg-Landau (TDGL) [3] equa- 
tions on the dynamic of the OP in the PSG has been 
established [j-lS|- Besides, it has been shown that in 
the limit of strong pair breaking effect due to interaction 
with phonons the GVG demonstrates strong hysteresis in 
the constant current regime [^ . All these results (see also 
[7|,|8[) demonstrate that different time-periodic and time- 
quasiperiodic solutions arise with the change of the volt- 
age or current [9] . Recent experimental observations have 
revealed the space-time arrangement of the PSGs. Using 
the low temperature laser scanning microscopy technique 
it was demonstrated that each voltage jump on the CVC 
corresponds to generation of a new PSG [10]. The au- 
thors have observed the creation of one PSG at a certain 
critical current. Further increase of the current leads 
to the spatial rearrangement of the PSCs - two PSCs 
appear symmetrically with respect to the center of the 
channel. Finally, the third PSG appears in the middle of 
the channel (see Fig. 2 of Ref.[10]). All these changes of 
the PSCs' arrangement immediately induce the change 
of the CVC's slope. However, some properties of the re- 
sistive state of a superconductor are not investigated. 

Our paper is structured as follows. First, we introduce 



the TDGL equations. Then we estimate analytically the 
averaged voltage and the period of the oscillating solu- 
tion for the relatively small currents. Finally, we analyze 
the bifurcation points of the TDGL equations and reveal 
their properties. 

The TDGL model is described by the equations: 
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Here ^p = p exp{iO) is the complex OP, where p and 
are the modulus and phase of the OP respectively. The 
distance and time are scaled in units of the coherence 
length ^ and phase relaxation time tq = A-kX^cFu/c^ re- 
spectively, where A is the penetration depth, Gn is the 
normal state conductivity, and c is the speed of light. The 
electrostatic potential (j) is measured in units of 0o/2'7rcr^, 
where ^o = ^hc/e is the flux quantum, e is the electric 
charge, h is the reduced Planck constant. The current 
density j is defined in units of (^oc/Stt'^X'^^. The only 
parameter left is u = Tp/ro^ where Tp is the relaxation 
time of the amplitude of the OP. From the other side 
u = i^ /l%^ where Ie is the penetration depth of the elec- 
tric field. The penetration depth increases if the inelastic 
scattering time on phonons is large enough. It can be ap- 
proximated in the case of gapless superconductivity if we 
assume that the phenomenological parameter u is small 
u <l 
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In the present paper we are considering the supercon- 
ducting wire of a length Lq/^ = 10.88 with the follow- 
ing boundary conditions: p{—L/2) = p{L/2) = 1 and 
d(j){-L/2)/dx = d(t){L/2)/dx = 0. The absence of the 
electric field at the end of the wire determines the gradi- 
ent of the phase: dO{-L/2)/dx = dO{L/2)/dx = j. 

Let us start with the features common to the cases 
of the channels of lengths L = 2Lo,4Lo, studied in Ref. 
[12]. The second TDGL equation in the steady state 
j = k{l — k'^) possesses two roots when j < jc- One 
of them is stable and the other is unstable. At j = jc 
they collide with each other. It is accompanied by the 
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FIG. 1: The period of the solution as a function of the current. 
The sohd hne represents the result of Eq. (3). Arrow indi- 
cates the period-doubling bifurcation point jd. Insets repre- 
sent projection of the limit cycle trajectory to the (yo(0), ^(0)) 
plane before (a) and after (b) the bifurcation point. Here p(0) 
is the modulus of the OP, ^(0) is the electric field both in the 
center of the wire. 



creation of a limit cycle. Besides, only one Lyapunov 
exponent crosses zero at this moment [6]. Therefore, the 
first singularity of the CVC at j = jc is a saddle-node 
homoclinic bifurcation p^. The period of oscillations 
in the vicinity of this bifurcation is determined by the 
formula 11 2l: 
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We have solved Eqs.(l) and (2) numerically with the 
help of the fourth order Runge-Kutta method. Tak- 
ing into account the finiteness of the channel, we have 
obtained that the critical current in our system jc is 
not equal to the one in the Ginzburg-Landau theory 
JGL = 2/3 a/3 for an infinite wire. Our calculations show 
that the critical current jc/JGL = 1.016. 

The results of our calculations for the period of the 
oscillating solution together with the analytical estima- 
tion (Eq.(3)) are presented in Fig.l. There is a very good 
agreement of numerical results with Eq.(3) over two or- 
ders of magnitude in (j — jc)/ jc 

From the Josephson relation [l| and Eq.(3) it follows 
that the voltage in the vicinity of the saddle- node homo- 
clinic bifurcation is proportional to the square root of the 
current V ex [(j— jc)/jc]^^^- Such behavior of the voltage 
for the region of currents jc < j < jd is represented on 
Fig. 2d. For this range of currents the PSC appears in 
the middle of the channel (Fig. 2a) in agreement with 
the experimental results [IQj. 

Further increase of the current causes the bifurcation of 
the periodic solution. It causes an increase of the period 
of the limit cycle at j = jd. On the insets to Fig.l we 
have plotted the phase portraits of the system before (in- 
set a) and after (inset b) the bifurcation point in coordi- 
nates {p{x = 0), E{x = 0)). A single-loop (period-1) limit 



FIG. 2: The CVC of the channel. Insets a-c show the space- 
time arrangement of the PSCs in the different regions of the 
CVC. Inset d represents the CVC for j < jd in comparison 
with analytical formula (see the text). 



cycle transforms to a double-loop (period-2) limit cycle 
at jci/JGL = 1.125. As it follows from Fig. 2b, the space- 
time arrangement of the PSCs is changed: two adjacent 
PSCs are shifted in the opposite directions with respect 
to the center of the wire in agreement with the experi- 
mental observations [lOj. Now the period includes two 
PSCs. Therefore, the type of the singularity which ap- 
pears at j = jci is the period-doubling bifurcation. As a 
result of the period-doubling, a new frequency UJ2 = cJi/2 
appears in the spectrum of an electromagnetic radiation 
generated by the current. It should be noted that the 
presence of the period-doubling bifurcation in the su- 
perconducting channels in the voltage driven regime has 
been reported earlier [9|. However, the authors observed 
the PSCs in the center of the channel, in contrast to our 
results. 

The last, but not the least bifurcation is the destruc- 
tion of the limit cycle. At jc2/ JOL = 1.44 the limit cycle 
looses stability. In this case two PSCs situated symmet- 
rically with respect to the center of the wire are accom- 
panied by the third one appearing in the center of the 
channel (see Fig. 2c). Such arrangement of the PSCs is 
confirmed by the experimental results [lO|. This situa- 
tion is in contrast to the longest channels cases, where 
we have obtained the appearance of the third PSC after 
destruction of the chaos (L = 4I/o) or inside the range 
of currents corresponding to the chaotic behavior of the 
system (L = 2Lo) [12]. Here we observe the appearance 
of the chaos and the third PSC simultaneously. We have 
observed a voltage discontinuity at this bifurcation point. 
For certain currents from the area j > jc2 we have ob- 
served so-called "periodic windows", areas of the param- 
eter values, where the limit cycle with a relatively large 
period becomes stable [14]. However, it does not induce 
any singularity on the CVC. For the range of currents 
j > jc2 the solution remains oscillating, but it looses 
the periodical properties. The Fourier transform of the 
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the spectrum. Finally, the appearance of the third PSC 
causes the destruction of the limit cycle. 
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FIG. 3: The Fourier amplitude of the voltage as a function of 
the frequency for a) j < jc2, b) j > jc2. 



voltage for j > jc2 possesses a lot of different frequen- 
cies (see Fig. 3b) instead of J-function like spectrum (see 
Fig. 3a) before the singularity. The same behavior can 
be observed for the cases of the longer channels, studied 
inRef. ^. 

To conclude, we have explored the properties of the 
narrow superconducting channel. We have demonstrated 
that the singularities of the C VC correspond to a number 
of different bifurcation points of the TDGL equations. 
The voltage appearance is accompanied by the saddle- 
node homoclinic bifurcation. It causes the formation of 
the limit cycle with a diverging period when j -^ jc. 
The voltage V ex {j— 2c)^l^ in this region. We have also 
analytically estimated the period of oscillations in the 
vicinity of this bifurcation point. The second singular- 
ity corresponds to the period-doubling bifurcation. As a 
result of this bifurcation, a new frequency equals to the 
half of the frequency before the bifurcation appears in 
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